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1.  Introduction 

In  [2]  conditions  were  given  under  which  it  could  be  proved  that  the 

problem 

inf  f(x) 
x  e  C  C  Rn  , 

where  C  is  a  convex  compact  set,  had  only  one  local  infimum  which  was 
taken  on  at  a  unique  point  x*  in  C  .  The  solution  was  shown  to  be  uncon¬ 
strained  in  the  sense  that  the  1  by  n  derivative  vanished  there,  i.e., 

f'(x*)  =  0  . 

Furthermore,  the  amount  by  which  the  function  f  evaluated  at  some  point 
in  C  exceeded  its  global  minimum  value  was  given  by  an  explicit  formula. 

In  this  paper,  similar  results  are  obtained  for  the  equality  constrained 
and  then  the  inequality  constrained  problem.  In  Section  2  some  preliminary 
lemmas  are  developed  for  the  problem  of  solving  simultaneous  nonlinear  equa¬ 
tions.  In  Sections  3  and  4  bounds  are  developed  for  an  isolated  local  con¬ 
strained  minimizer  in  C  ,  and  on  the  minimum  value  associated  with  that 
minimizer.  It  is  also  shown  that  if  slightly  stronger  conditions  are  placed 
on  the  problem  (verifiable)  then  that  local  minimizer  is  a  global  one  in  C  . 


I 
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2.  Solving  Nonlinear  Equations — Preliminary  Lemmas 

Let  a(z)  be  an  Nxl  vector  function  which  is  continuously  dif- 

N 

ferentiable  in  an  open  set  containing  a  given  compact  convex  set  C  CR 
Assume  that  a'  (z)  ,  the  NxN  derivative  matrix  of  a(z)  ,  has  an  in¬ 
verse  at  every  point  in  C  .  Let  denote  the  interior  of  C  .  De¬ 
fine  U  =  {u|u=a(z)  for  some  z  e  C}  .  Consider  any  point  z  e  .  Put 

u  =  a(z)  .  It  follows  from  the  inverse  function  theorem  that  there  ex- 

/V  0  /\  ^  A 

ists  an  open  set  C  C  C  containing  z  ,  an  open  set  U  C  U  containing  u  , 

and  a  continuously  differentiable  function  g(*)  defined  on  U  such  that 

g[a(z)]  =  z  ,  for  z  in  6  ;  a[g(u)]  =  u  ,  for  u  in  U  . 


Lemma  1.  Suppose  for  some  z^  e  C  , 

NT  /  \  = 


Then: 


N(Zq)  =  (Zq  -  a'  (y)  a(zQ)  |  yeC}  C  C 


a(zg)t  e  U  ,  for  all  t  such  that  O^t^l  , 


and  there  exists  a  point  z*  in  N(z^)  such  that 


a(z*)  =  0  . 


Proof:  Let  t  be  the  smallest  value  greater  than  or  equal  to  zero  such 
that  g[a(Zg)t]  is  defined  for  O^^t^l  .  (Clearly  t<l  ,  since 

Zq e  C^  and  because  of  the  inverse  function  theorem.)  Assume  t  is  such 


that  0  <  t  <  t  <  1  .  Now 


dg[a(z) ]/da  -  a' (z)  1  . 


Thus, 


g[a(zQ)t]  *  g[a(zQ)]  +  /  (dg[a(z0)(l-s) + (zQ)st] /da}ds  a(zQ)(t-l) 

1  _l 

■  z 0  -  /  a' [g{a(zQ) (1-s+st) }]  ds  a(zQ)(l-t)  . 


(Note  because  0  <  t  <  1-s+st  <  1  ,  g{a(Zg) (1-s+st) }  is  defined  and  in  C  .) 


■> 


Now, 
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y  =  zQ  -  Iq®'  tg{a(zQ)(l-s+st)}]  ^s  a(zQ) 


(5) 


is  in  C  .  To  see  this,  assume  the  contrary.  Because  C  is  a  convex  set, 

A  A  A 

there  would  exist  a  hyperplane  (a, 3)  which  separates  y  from  C  with  the 
properties  that 


and 


az  -  3  £  0  for  all  z  e  C  , 

if  az  -  3  >  0  z  i  C  , 

ay  -  3  >  0  . 


(6) 


Using  (5)  and  the  mean  value  theorem,  there  is  an  s  where 
0  <  s  <  1  such  that 


a[zQ-a' [g{a(z0)(l-s+st)}]-1a(z0)]  -  3  >  0  . 

_  A  A  A  A 

(Note  because  t<l-s+st£l  ,  g{a(zQ) (1-s+st) }  is  defined  and  in  C  .) 
Using  (6)  this  contradicts  (1).  Thus  y  e  C  .  Now  zQt  +  y(l-t)  gives 
(4) .  Since  C  is  a  convex  set  and  0  <  t  < 1  ,  and  zQ  e  C°  ,  (4)  is  also  con¬ 
tained  in  C  .  Since  a' (z)  is  continuous  it  follows  that  g[a(Zg)t]  e  C  . 
If  t  >  0  the  same  reasoning  as  above  implies  g[a(zQ)t]  e  C°  ,  and  the 

inverse  function  theorem  implies  the  existence  of  values  less  than  7  and 
close  to  it  for  which  a(zQ)t  e  U  .  By  contradiction  then,  7=0.  This 

proves  (2)  and  (3). 
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3.  The  Equality  Constrained  Problem 

The  equality  constrained  nonlinear  programming  problem  can  be  written 


minimize  f(x) 


subject  to  x  e  R  =  tx|h.(x)  =  0  ,  for  j=l,...,p}  and  x  e  C 

J  ^ 


(7) 


where  C  is  a  compact  convex  subset  of  R 
x 


n 


The  functions  f  ,  {h^}  are 


assumed  to  be  twice  continuously  differentiable  in  an  open  set  containing  . 

The  vector  h(x)  is  the  p  by  1  vector  whose  jth  component  is  h^ (x)  .  The 

following  development  will  be  general  in  that  the  particular  numerical  method 
for  generating  the  null  space  of  h'(x)  (the  p  by  n  derivative  matrix 
of  h(x))  will  not  be  specified. 


Suppose  S(x)  denotes  a  matrix  function  which  gives  the  null  space 
of  h' (x)  .  That  is,  a  necessary  and  sufficient  condition  that  h'(x)z  =  0 
is  that  z  =  S(x)v  for  some  v  .  Also  associated  with  h' (x)  is  a 
pseudo  inverse  h'(x)//  ,  i.e.,  a  matrix  satisfying 


h'(x)  h'(x)#h'(x)  =  h'(x)  . 

Note  that  when  the  rank  of  h’ (x)  is  p  ,  h’(x)h’ (x)#  =  I  . 

It  is  assumed  only  that  S(x)  and  h'(x)//  are  related  in  that 
[I-h' (x)//h' (x) ]  =  S(x)W(x)  for  some  matrix  W(x)  . 

The  Lagrangian  function  associated  with  (7)  is 

L(x,a)  =  f(x)  -  aTh(x)  . 

T  T  T 

The  vector  z  of  Section  2  is  to  be  identified  with  z  =  (x  ,a  )  , 


is  (N»(n+p))  by  1.  The  following  result  is  easily  obtained: 


5 


a'(z)  = 


»°)  >  -h’(x)T\ 

-h '  (x)  ,  0  /  • 


Let  be  the  convex  hull  of  {f 1  (x)h'  (x) //  |x  e  C^}  .  Define 

C  =  CxX  Ca  .  Define  H(x,a)  =  S(x)VxL(x,a)S(x)  . 

Theorem  1.  Suppose  for  every  z  e  C  , 

(i)  h'(x)  has  full  rank  (equal  to  p) 

(ii)  S(x)  is  a  continuously  differentiable  n  by  (n-p) 

matrix  function  generating  the  null  space  of  h'(x)  , 
and  that 

(iii)  H(x,a)  is  a  positive  definite  matrix. 

Then  a' (z)  ^  exists  for  all  z  e  C  and  is  given  by 


a'  (z) 


rs(x 

-1  = 

-h' 


S(x)H(x,a)  1S(x)T,-[I-S(x)H(x,a)  (x)TVxxL(x,a) ]h' (x)#“ 


(x)#T[I-VxxL(x,a)S(x)H(x,a)  1S(x)T] ,F 


where 


F  =  -h'(x)#T[V*L(x,a)-V*  L(x,a)S(x)H(x,a)_1S(x)V  L(x,o)  )h' (x)#  . 

XX  XX  XX 


Let  g(*)  denote  the  inverse  function  as  defined  in  Section  2.  The 
notation  S(t)  will  mean  S[g{a(zQ)t}]  and  similar  notation  used  for  the 

other  quantities.  Suppose  xQ  e  Cx  is  given.  Define  aQ  =  [f '  (x^h'  (xQ)//]T  , 

and  set  z*  =>  (x^.a^)  .  Define  N(zQ)  =  {zQ-a' (y)-1a(z0) | y  e  C}  .  Assume 

further  (iv)  that  N(zQ)  C  C  and  that  zQ  e  C°  . 

Then  there  exists  a  z*  e  N(zQ)  such  that  ?xL(x*,a*)  =  0  ,  h(x*)  =  0 
is  an  isolated  local  minimizer  to  Problem  (7),  and 


(10) 


f(x*)  -  f(xQ)  =  /5{tVxL(x0,a0)TS(t)H(t)“1S(t)TVxL(x0,a0) 

+  f’(t)[I-S(t)H(t)_1S(t)TVxxL(t)]h,(t)#h(x0)}dt  . 

Proof :  Because  H(x,a)  is  positive  definite,  its  inverse  exists.  The 
remainder  of  the  first  assertion  above  can  be  verified  by  multiplying 
the  two  matrices  together. 

The  existence  of  z*  is  guaranteed  by  Lemma  1.  Since  VxL(x*,a*)  =  0  , 

(x*,a*)  satisfies  the  first  order  necessary  conditions  for  an  equality  con¬ 
strained  extremum.  Assumption  (iii)  means  then  that  the  second  order  sufficiency 
conditions  for  an  isolated  local  minimizer  are  satisfied  at  x*  . 

From  Lemma  1  it  also  follows  that  g{a(zg)t}  is  defined  and  in  C 
for  all  0  <_  t  <_  1  .  Define  the  composite  function 

F(t)  =  f[g{a(zQ)t}]  . 

Then 

f(x*)  =  F(0)  =  F(l)  +  /°  F'(t)dt,  where  F(l)  =  f(xQ)  . 

Using  the  chain  rule  of  differentiation, 

F’(t)  =  df [g{a(zQ)t}]/dt  =  df [g{a(zQ)t}]/dz  •  dg{a(zQ)t}/da  •  d[a(z0>t]/dt 

=  f  ’  (t)S(t)H(t)"1S(t)TVxL(x0,a0)+f ’ (t)  [I-S(t)H(t)-1S(t)  VxL(t)  ]  - 

h’(t)//h(x0)  . 


Now  h’(t)S(t)  *  0  ,  so 


f'(t)S(t) 


[f'(t)-a(t)h'(t)]S(t)  =  VxL(t)TS(t) 


tVxL(x0,a0)TS(t)  .  (11) 


Using  this  fact  completes  the  proof  of  the  theorem. 


4.  The  Inequality  Constrained  Problem 


The  notation  of  the  previous  section  will  be  used  and  the  preceding 
results  extended  to  the  inequality  constrained  problem 

minimize  f(x) 

subject  to  x  e  R  =  {x  |  lv.(x)  >;0  ,  for  j=l,...,m)  ,  (12) 

and  xe  C  where  C  is  a  compact  convex  subset  of  Rn  . 

xx 

Lemma  2.  If  x*  is  an  isolated  local  minimizer  for  the  Problem 
minimize  f(x) 

subject  to  xe  {x|h.(x)=>0  ,  for  j=l, .  • .  ,p^m}  H  C  ,  (13) 

3  % 

and  if  x*eR  ,  then  x*  is  an  isolated  local  minimizer  for  Problem  (12). 
(Note:  the  proof  is  obvious  and  therefore  omitted.) 

In  what  follows,  conditions  will  be  given  which  guarantee  the 
existence  of  an  isolated  local  constrained  minimizer  in  C  for  the 
Problem  (13).  Without  loss  of  generality  it  will  be  assumed  that  the 
vector  h(x)  will  consist  of  the  first  p  functions. 

Theorem  2.  For  the  inequality  constrained  Problem  (13)  assume  as  in 
Theorem  1.  Then  there  exists  a  z*  e  N(z q)  such  that  v  L(x*,a*)  =  0  , 

and  h(x*)  =  0  . 

Proof :  The  proof  of  this  theorem  parallels  that  for  Theorem  1  and  will  not 
be  given.  Note  that  it  is  not  possible  to  conclude  that  x*  is  a  local 
minimizer  for  the  inequality  constrained  problem.  Several  more  conditions 
must  be  placed  on  the  problem  quantities. 

Corollary  1.  Assume  as  in  Theorem  2.  If  at  x*  ,  a*  >  0  and  x*  e  R  , 
then  x*  is  an  isolated  local  minimizer  for  the  Problem  (12) . 

Proof:  The  proof  follows  because  under  the  assumptions,  x*  satisfies  the 
second  order  sufficiency  conditions  for  a  constrained  minimizer  to  (13). 
Using  Lemma  2,  the  corollary  is  proved. 

5.  Conditions  Under  Which  the  Isolated  Local  Minimizer 
Is  a  Global  Minimizer  in  the  Convex  Compact  Set 

Some  notation  needs  to  be  established.  Let  x  be  any  point  in 
Cx  n  R  .  Define  for  0  <_  t  <_  1  ,  y(t)  ■  x*(l~t)+xt  ,  and 

y(a,t)  ■  x*(l-s)+[x*(l-t)+xt]s  ,  for  0  <_  s  <_  1  . 
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Theorem  3.  Assume  as  in  Theorem  1.  In  addition  the  following  assumptions 


are  made. 


(i)  /J  h'[y(t)]dt  has  rank  p  , 

(ii)  vT/q  /q  V^xL[y(s,t),a*]dsdt  v  >  0  for  all 
v  ^  0  where  /Jh' [y(t)]dt  v  =  0  , 


/ o  /J  V^xL[y(s,t),a*]dsdt  ,  -  jJ  h' [y  (t)  ]Tdt' 
“  /q  h' [y(t)]dt  ,  0  , 


has  an  inverse.  Denote  this  by 


(:::) 


T  T 

\  where  A  is  n  by  n  ,  B  is  n  by 


r  1  T 

If  (iii) »  “  Jq  f ' [y(t) ]  dt  B  >  0  ,  (component  by  component)  then 

f (x)  >  f(x*)  .  If  (i),  (ii),  and  (iii)  hold  for  all  x  e  C  ,  then  x*  is 

x 

the  unique  global  minimizer  for  £  in  Rfl  C 

x 


Proof:  Consider  the  following  (recall  VxL(x*,a*)  =  0  ,  h(x*)  =  0): 

// J  VxL[y(t),a*]\  /0\  //J  /J  V^xL[y(s,t),a*]dsdt  ,  -  /J  h' [y (t) ]Tdt\  /Ax 


-  h'(x) 


-  /J  h'[y(t)]dt  , 


where  Ax  =  (x-x*)  ,  Aa  =  (a*-a*)  =  0  . 


The  existence  of  the  above  inverse  is  guaranteed  by  the  theorem 
hypotheses.  It  can  be  written  explicitly  after  the  following  notational 
considerations  are  given.  Let  S  be  any  n  by  (n-p)  matrix  generating  the 

null  space  of  N  =  /^  h'[y(t)]dt  .  Let  N//  be  a  psueudo  inverse  &•£ 

related  to  S  in  that  there  exists  a  matrix  W  such  that  [I-N//N]  =  SW  . 
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1.  X  A 

Let  H  -  /q  f o  ^xxL[y(s,t) ,a*]dsdt  .  The  explicit  inverse  (which  is 

invariant  to  the  choice  of  the  particular  null  space  matrix  S  and  pseudo 
inverse  N//)  is  then 

/  H,  -NT\_1  1  S(STHS)-1ST  ,  -[I-S(STHS)'1STH]N// 

\-N,  0  /  |-N#T[I-HS(STHS)-1ST]  ,  -  n#t[h-hs(sths)_1sth]n* 

Now 

f(x)  -  f(x*)  =  /£  f'[y(t)]dt  Ax 

=  /£  f ’ [y (t) ]dt  S(STHS)-1ST/J  VxL[y(t),a*]dt 

-  /£  f’[y(t)]dt  BTh(x) 

=  /£  VxL[y(t),a*]TdtS(STHS)_1ST/£  VxL[y(t),a*]dt 

-  /£  f’ [y(t) ]dtBTh(x)  . 

Since  x  e  R  ,  h(x)  0  .  Thus  both  terms  on  the  right  above  are 
X  X 

nonnegative.  If  ^xL[y(t),a*]  dtS^O  ,  the  first  term  above  is  >  0  . 

The  second  term  is  >0  if  h^ (x)  >  0  for  one  of  j  =  l,...,p  .  The 

X  X 

only  way  both  terms  can  be  equal  to  zero  is  if  /  VxL[y(t) ,  a*]  atS  =  0 
and  h(x)  =  0  . 

If  both  these  conditions  hold  it  is  easy  to  show  that  Ax  =  x  -  x*  =  0 
Thus  f(x)  >  f(x*)  when  x  t  x*  . 

Several  interesting  results  come  out  of  this  analysis.  There  are 
two  estimates  of  the  generalized  Lagrange  multipliers  associated  with  a 
local  minimizer.  The  first  one  is 


“o  =  f,(xo)h'(xo^  * 


(15) 
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This  is  used  in  the  computation  of  the  Lagrangian  Hessian  matrix.  The 
second  estimate  is 

f' (xQ)  [I-S(x0)H(x0,a0)  1S(x0)TV^xL(x0,a0]h' (xQ)//  .  (16) 

This  is  the  one  which  determines  whether  or  not  a  constraint  is  to  be 
considered  binding  at  the  local  solution.  The  latter  estimate  is  really 
the  estimate  of  the  multipliers  (15)  at  the  minimizer  of  f  in  the  sub¬ 
space  h(x)  =  Mxq)  •  In  algorithm  development,  then,  the  choice  of 

whether  or  not  to  include  a  constraint  in  the  set  of  those  considered 
active  should  be  made  on  the  basis  of  (16)  rather  than  (15) . 


Discussion  and  Example 


The  importance  of  the  second  order  estimate  of  the  multipliers  was 
first  pointed  out  by  Gill  and  Murray  in  their  seminal  paper  [1] .  The 
estimate  in  (16)  above  corresponds  exactly  to  their  "refinement"  given 
in  (7.4)  of  [1].  Their  vector  p  there  is  the  approximation  required 
to  find  the  minimizer  of  the  function  in  the  equality  constrained  sub¬ 
space.  In  the  current  paper  this  is  the  transpose  of  the  vector 

-f'(x0)S(x0)H(x0,a0)  ^(Xq)7 

in  (16)  above. 


The  quantity  (10)  which  gives  the  exact  difference  between  the 
value  of  the  objective  function  at  the  minimizer  in  the  space  where 
h(x)  =  0  and  its  value  at  some  point  x^  may  be  considered  (for  intui¬ 
tive  explanatory  purposes)  to  consist  of  two  parts.  Let  x  denote  the 
minimizer  for  the  problem:  minimize  f(x)  subject  to  h(x)  *  h(x^)  » 

xE  .  The  first  term  under  the  integral  in  (10)  is  an  approximation 

to  f(x)  -  f(Xg)  •  It  is  well  known  that  the  difference  in  the  optimal 

objective  function  value  resulting  from  a  perturbation  in  the  right-hand 
side  of  the  constraints  is  approximated  by  the  sum  of  the  Lagrange  multi¬ 
pliers  times  the  perturbations.  The  vector  premultiplying  h(xg)  is  the 


* 
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integral  of  a  vector  of  second  order  multiplier  estimates  at  unconstrained 
minimizers  of  the  form  (16).  This  can  be  regarded  as  the  appropriate 
vector  of  Lagrange  multipliers  needed  to  approximate  f(x*)  -  f(x)  . 

Thus  (10)  says  (in  effect)  that  f(x*)  -  f(x^)  =  [f(x)-f(xQ)]  + 

[f(x*)-f(x)]  . 

In  terms  of  the  inverse  function  approach  x  consists  of  the 
first  n  components  of 


The  exact  difference  f(x^)  -  f(x)  can  be  computed  by  integrating  from 

(xJ,aJ)T  to  (xT,aT)T  along  g[ (tVxL(x0>a0)T,  h(xQ)T)T]  .  The  exact 

second  term,  f(x*)  -  f(x)  ,  can  then  be  obtained  by  integrating  from 

(xT,aT)T  to  (x*T,a*T)T  along  the  curve  g[(0^,  h(x^)T)T]  .  The  reason 

that  the  current  approach  is  taken  is  that  the  existence  of  the  inverse 
function  along  the  "diagonal  path"  can  be  checked  by  (1) .  The  piecewise 

two-path  segment  described  above  may  not  exist  in  the  region  C  x c 

x  ct 

without  stronger  assumptions. 

The  practical  application  of  this  theory  is  that  from  (1),  it  is 
possible  to  place  upper  and  lower  bounds  on  each  component  of  the  mini- 
mizer  and  multipliers,  and  from  (10)  it  is  possible  to  place  an  upper  and 
lower  bound  on  the  optimal  objective  function  value.  It  is  in  general 
not  possible  to  compute  the  exact  set  N(Zq)  of  (1)  or  the  right-hand 

side  of  (10).  Using  numerical  techniques,  one  can  compute  upper  and 
lower  bounds  on  the  desired  quantities .  For  the  purely  unconstrained 
case,  the  computational  requirements  were  discussed  in  [3].  The  con¬ 
strained  case  is  more  complicated  and  a  future  paper  will  discuss  the 
linear  algebra  required  to  do  this. 

Essentially  what  has  to  be  done  is  to  replace  each  scalar  function 
required  in  the  formulas  by  two  numbers,  an  upper  and  a  lower  bound  on 
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the  range  that  that  scalar  function  can  take  when  x  is  restricted  to 
,  and  a  to  C  .  Techniques  required  to  do  this  have,  to  a  cer¬ 
tain  extent,  been  developed  into  a  general  theory  of  "interval  arith¬ 
metic"  [4],  Every  element  of  a  vector  and  matrix  can  be  systematically 
replaced  by  two  numbers.  The  problem  becomes  more  difficult  when  an 
interval  inverse  matrix  is  required,  i.e.,  a  matrix  of  interval  ele¬ 
ments  which  contains  the  inverse  of  all  possible  matrices  implied  by  the 
interval  matrix.  In  order  to  get  a  set  bounding  (1) ,  an  interval  inverse 

matrix  is  required  for  a' (y)  *  ,  where  the  form  taken  for  the  constrained 
problem  is  given  by  (8).  In  the  following  example,  ad  hoc  procedures  are 
used  to  get  the  interval  bounds.  A  few  of  the  details  are  given,  but  the 
bulk  of  the  computations  must  be  taken  on  faith. 

The  rules  for  computing  intervals  for  the  range  of  functions  (or 
interval  extensions  as  they  are  sometimes  called)  are  very  simple.  The 
notation  [a,b]  is  used  to  indicate  that  an  upper  bound  on  the  range  of 
a  scalar  function  is  b  and  a  lower  bound  is  a  .  If  two  scalar  valued 
functions  have  interval  extensions  [a,b]  and  [c,d]  ,  respectively, 
then  the  sum  has  the  bounds  [a+b,c+d]  .  The  range  on  the  product  is 
[min(ac,ad,bc,bd) ,  max[ac,ad,bc,bd) ]  .  Subtraction  is  similar  to  addi¬ 
tion.  Division  is  similar  to  multiplication  as  long  as  the  divisor  does 
not  contain  zero.  The  other  major  requirement  is  that  interval  extensions 
be  computable  for  functions  of  a  single  variable  operating  on  an  interval. 
In  one  sense  this  can  be  viewed  as  an  optimization  problem.  For  most 
functions  of  a  single  variable,  the  min  and  max  of  that  function  over  an 
interval  are  easily  discernible.  Suppose  a  scalar  t  is  known  to  have 
a  range  of  [a,b]  .  Then  the  scalar  function  t^  has  a  range  of 

[min(a^,b^),  max(a^,b^)]  if  0  i  Ia,b]  ,  and  [0,  max(a^,b^)]  if 
a^0<b  .  A  scalar  function  which  has  the  range  [a, a]  will  in  the 
following  be  indicated  by  a  for  simplicity.  Thus,  the  notation  a[c,d] 
will  indicate  the  range  of  the  product  of  two  functions,  one  of  which  has 
range  [a, a]  ,  the  other  [c,d]  .  The  resulting  range  is  [min(ac,ad), 
max(ac,ad)]  .  The  major  research  area  needing  to  be  investigated  is 
obtaining  an  Interval  extension  of  a  matrix  inverse. 
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The  problem  is 

minimize  -x  +  y 

2  2 

subject  to  h  =  -x  -y  +  1  >_  0  , 

C1  =  { (x,y)  |*2  £  x  <  2  ,  -  2  <  y  <  -  .2}  . 


The  quantities  required  are 

h'  =  (-2x,  -2y)  ,  (h')l 


•n 


,  2.-1  n2  _ 

(x  +y  )  ,  V^L 


0 


[2(x2+y2) ]_1 


2  2 

The  multiplier  estimate  is  given  by  f'(h’)#  *  (x-y)/[2(x  +y  )]  .  The  top 

2  2 

component  of  (h')#  is  -x/[2(x  +y  ) ]  .  One  can  show  that  the  minimum  of 
this  function  in  the  rectangle  is  taken  on  at  x  =  .2,  y  ■  -.2  and  the 
maximum  at  x  =  2,  y=~2.  Thus  the  range  of  the  first  component  is 

[-1.25,  -.125]  .  Similar  considerations  yield  the  range  of  the  second 
component  to  be  [.125,  1.25]  .  The  range  on  a  is  therefore 
(-1)1-1.25,  -.125]  +  (1)[.125,  1.25]  ■  [.25,  2.5]  .  This  then  is  the  set 

Ca  ‘ 


The  matrix  whose  inverse  is  desired  [from  (8)]  is 


The  interval  extension  of  this  inverse  when  x,y  ,  and  a  are  restricted 
as  indicated  above  is 
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B 


([.001980,  1.98020] 
[.019802,  1.0] 
[.024752,  1.25] 


[.019802,  1.0] 
[.001980,  1.98020] 
[-1.25,  -.024752] 


[.024752,  1.25]  \  (17) 

[-1.25,  -.024752]  j  . 
[-15.625,  -.015625]/ 


Let  (x0,y0)  ■  [-7,  -.71]  .  Thus  aQ  =  .709184,  hQ  =  .0059,  and 

=  [-.007142,  -.007042]^  .  Then  N(Zq)  of  Theorem  2  is  contained  in 


([  .700300,  . 728560]\  /[  0.2,  2]  \ 

[-.717220,  -.689060]  J  C  (  [-2  ,  -0.2]  J  . 

[  .608371,  .717845]  /  \[  0.25,  2.5]/ 


which  is  clearly  in  the  original  rectangle  of  bounds. 


The  bound  on  the  optimal  objective  function  value  is  calculated 

using  (10).  The  first  term  under  the  integral  involves  S(t)H(t)  ^S(t)  . 
Note  from  (9)  that  this  occurs,  for  any  (x,a)  ,  in  the  upper  left-hand 
corner  of  the  inverse  of  a'  .  This  is  thus  available  from  (17)  above. 

An  interval  extension  (for  every  t)  of  V  L(x_,a~)TS(t)H(t)  *S(t)TV  L(xrt,art) 

X  0  U  x  0  U 

is 


(-.007142,  -.007042) 


([.001980,  1.98020], 
[0.19802,  1.0]  , 


[0.19802,  1.0] 
[.001980,  1.98020] 


X-. 007142 
-.007042 


=  [.000002191,  .00029976]  . 


A  term  of  one-half  comes  from  the  integration  of  a  constant  times  t  for 
t  between  0  and  1.  Thus  the  contribution  of  the  first  term  is 
[.0000010955,  .00014988]  . 

The  second  term  under  the  integral  in  (10)  involves  the  negative 
of  the  matrix  which  occurs  in  the  upper  right-hand  corner  of  (9) .  This 
needs  to  be  premultiplied  by  the  objective  function  derivative  (-1,  1) 
and  postmultlplied  by  Mxq)  to  yield 
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(-[  .024752,  1.25]  \ 

)  (.0059)  =  [.00029208,  .01475] 

-[-1.25  ,  -.024752]/ 

Adding  the  two  contributions  yields  the  bounds 
[.0002931755,  .01489988]  . 

Since  -  —1.41  ,  the  value  of  f*  must  lie  in  the  interval 
[-1.42490,  -1.41029]  . 
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